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We show that a wide range of spin clusters with antiferromagnetic intracluster exchange inter- 
action allows one to define a qubit. For these spin cluster qubits, initialization, quantum gate 
operation, and readout are possible using the same techniques as for single spins. Quantum gate 
operation for the spin cluster qubit does not require control over the intracluster exchange interac- 
tion. Electric and magnetic fields necessary to effect quantum gates need only be controlled on the 
length scale of the spin cluster rather than the scale for a single spin. Here, we calculate the energy 
gap separating the logical qubit states from the next excited state and the matrix elements which 
determine quantum gate operation times. We discuss spin cluster qubits formed by one- and two- 
dimensional arrays of s = 1/2 spins as well as clusters formed by spins s > 1/2. We illustrate the 
advantages of spin cluster qubits for various suggested implementations of spin qubits and analyze 
the scaling of decoherence time with spin cluster size. 



I. INTRODUCTION 

During the past years, the discovery of several pow- 
erful quantum algorithms^ has triggered substantial re- 
search efforts aiming at the implementation of a quantum 
computer in a physical system. The main difficulty is 
that qubits must be prepared, manipulated, and read out 
with high fidelity while decoherence is required to remain 
small. 2 Solid-state implementations of qubits exploit the 
versatility of nanoscale fabrication, but suffer from de- 
coherence times which are usually shorter than in many 
quantum optics proposals. 3 ElectrorA^ and nucleari^ 
spins have been identified as promising candidates for 
qubits in a solid state system. The main advantage of 
electron or nuclear spins is that they are natural two 
state systems and that decoherence times for the spin 
degree of freedorr&iS are usually larger than for charge 
degrees of freedom. 

Here we show that a wide variety of spin clusters are 
promising candidate systems for qubits. Qubits formed 
by several spins have so far mainly been discussed in 
the context of exchange-only quantum computing, iLiSii^ 
coherence-preserving qubits^, and quantum computing 
schemes in which the requirements on the control of ex- 
change interactions between spins are relaxed^ How- 
ever, all these schemes require control at the single-spin 
level, either with local magnetic fields^ or exchange in- 
teractionsiii*i^ For the spin clusters considered here, con- 
trol for both magnetic fields and exchange interactions is 
required only on the length scale of the spin cluster di- 
ameter. As we have shown in Ref. llfjl spin chains formed 
by an odd number of antiferromagnetically coupled spins 
s = 1/2 allow one to define a logical qubit. The logi- 
cal state of the qubit is encoded in the collective state 
of the spin cluster. Here we detail that this construction 
remains valid for a wide range of spin clusters, indepen- 
dent of the details of intracluster exchange interaction 
and spin placement. Initialization and readout of the 
spin cluster is achieved with the methods developed for 



single spins4£ The main advantage of spin clusters is 
that the requirements on spatial control can be traded 
for gate operation times. The scaling of the decoher- 
ence rate with the size of the spin cluster depends on the 
microscopic decoherence mechanism. While the decoher- 
ence rate induced by fluctuating local magnetic fields in- 
creases with cluster size, we show that magnetic dipolar 
interactions for the spin cluster qubit are smaller than for 
single spins. The optimum size of the spin cluster qubit 
is determined by a trade-off between the increase in gate 
operation times and the decoherence rate effected by lo- 
cal fluctuating magnetic fields, the decrease in magnetic 
dipolar interaction energy, and the relaxed conditions on 
local control. 

Any quantum computation can be decomposed into a 
sequence of one- and two-qubit quantum gatesiii For a 
single-spin qubit, the s z eigenstates | f) and | J.) are iden- 
tified as logical basis states |0) and |1), respectively^^ 
The phase shift gate can then be realized by a mag- 
netic field B z (t) and the one-qubit rotation gate C/ ro t by 
a transverse field B x (t) which rotates | f) into | J.) and 
vice versa. More generally, the equations 



\H'\0) = (l|#'|l)and (1\H'\0) ^0 



(1) 



constitute a sufficient condition that a Hamiltonian H' 
induces the unitary time evolution required for U TO % . For 
single spins, H' = g^BB x (t)s x fulfills Eq. QJ. Similarly, 
an exchange interaction = J„§i • §2 generates the uni- 
tary time evolution required for the square-root of SWAP 
gate- because, in the two-qubit product basis, 



(10|#»|01) ^ 0. 



(2) 



In contrast to a single spin s = 1/2, clusters formed by n c 
coupled spins are not intrinsically two-state systems. In 
order to prove that a logical qubit can be defined in terms 
of the energy eigenstates of a spin cluster we will (a) iden- 
tify spin clusters with a ground state doublet {|0), |1)} 
separated from the next excited state by an energy gap 
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A; (b) identify Hamiltonians H' and which satisfy 
Eqs. ||TJ and J2J and, hence, allow one to generate a uni- 
versal set of quantum gates; and (c) quantify leakage and 
decoherence for the spin cluster qubit. In particular, the 
evaluation of the matrix elements in Eqs. (JTJ and © and 
the quantification of excitation out of the computational 
basis (leakage) requires a detailed characterization of the 
states {|0), |1)} which is, in general, nontrivial. 

This paper is organized as follows. In Sec.^we discuss 
the computational basis states for spin-1/2 chains. For 
this simple geometry, it is possible to derive analytical ex- 
pressions for the matrix elements in Eqs. JQl and J5| for 
various anisotropics and spatially varying intracluster ex- 
change interaction. Section lTTll discusses the insensitivity 
of spin cluster qubits to the details of interactions within 
the cluster, such as the relative placement of spins and 
the exchange strengths. In Sec. II VI spins with spin quan- 
tum numbers larger than 1/2 are discussed. In Sec. Ivl 
we draw our conclusions. 



II. SPIN CHAINS 

For simplicity, we first consider a spin cluster qubit 
formed by a spin chain, 

n c -l 

H ^ ^ fj [J_L(sj,xSj-\-l,x & j ,yS j + ~h Jz^j^z^j-\-l,z] 

2 = 1 

where n c is odd and Jj_, J 2 > 0. The real numbers fj > 
account for a spatial variation of the exchange interac- 
tion, and J±fj (Jzfj) denotes the transverse (longitudi- 
nal) exchange interaction between sites j and j + 1. 



A. Isotropic spin chains 



For electron spins in quantum dots, the nearest neigh- 
bor exchange is usually of the Heisenberg form^ J = 
J± = J z . We first consider fj = 1, 



H = J 



(4) 



with J > 0. Note that this is an open spin chain; a closed 
spin chain would have a fourfold degenerate ground state 
multiplet for odd n c that would make it unsuitable for 
representing a single qubit. Because the intracluster ex- 
change interaction J is time independent and no exter- 
nal control is required, J can be adjusted already during 
sample growth. 

Spin chains have been studied in great detail during 
the past decades The theoretical description of 
the antiferromagnetic spin chain Eq. I0J is particularly 
challenging because the classical Neel ordered state is 
not an energy eigenstate and quantum fluctuations are 
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FIG. 1: Energy spectrum of an isotropic spin chain with n c = 
5 (left panel) and n c = 6 (right panel) . Energy eigenstates are 
sorted according to their quantum numbers of S z and their 
eigenenergies. 

pronounced. We define the operator of total spin, 

S a = ~y s,-. Q (5) 



for a — x,y,z. Energy eigenstates can be labeled ac- 
cording to their quantum numbers of total spin S and 
the ^-component of total spin, S z , because 



[H, S 2 ] = [H, S z ) = 0. 



(6) 



Due to the antiferromagnetic exchange coupling, states 
in which the total spin of the chain is minimized are en- 
ergetically most favorable^i For even n c , the minimum 
possible spin is S = 0, and the system has a nondegen- 
erate ground state. In contrast, for odd n c , there is a 
ground state doublet (Fig.QJi^i This parity effect is well 
known for thermodynamic quantities^ The energy gap 
A separating the ground state doublet from the next ex- 
cited state, 



" 2 mm 2n r 



(7) 



can be estimated from the lower bound of the des 
Cloiseaux-Pearson spectrum and the minimum wave vec- 
tor fcmjn = Tr/n c (Kei. E^) . Henceforth, we will restrict 
our attention to spin chains with odd n c . 

The requirements on a candidate system for qubits in- 
clude initialization of the quantum computer, a universal 
set of quantum gates, decoherence times long compared 
to gate operation times, and readout of the qubit. 2 



1. Definition of the spin cluster qubit 

For the chain with an odd number of sites [Fig. E{a)], 
we define the spin cluster qubit in terms of the 5* = 1/2 
ground state doublet by 



Sz\0) = -|0), 
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FIG. 2: Spin cluster qubit. (a) The energy eigenstates of 
an antiferromagnetic spin chain with an odd number of sites 
define the spin cluster qubit. (b) Wave function of |0) in 
the single-spin product basis for n c = 9. The size of each 
configuration is proportional to the probability of finding the 
corresponding product state in |0). (c) Similar to (b) but 
for state |1). (d) Spin density {ip\s itZ \ijj) for \ip) — |0) and 
\ip) — |1) (solid lines) for constant intrachain exchange cou- 
pling. Dashed lines indicate the corresponding data for spa- 
tially varying intrachain exchange with fj — sin(jn /n c ) . (e) 
Spin cluster qubits are coupled by a switchable interchain ex- 
change coupling J*(i). 



The states {|0), |1)} do not in general have a simple repre- 
sentation in the single-spin product basis, but rather are 
complicated superpositions of n c !/[(n c — 1 )/2] ! [(n c +l )/2] ! 
states [Figs. |2I[b) and (c)] as evidenced by the local mag- 
netization density [Fig. |2fd)] . The largest amplitude in 
this superposition corresponds to the Neel ordered states 
I T>i| l)a... I T)» c (|0)) and | |)r| T) 2 . . . | |)„ c (|1», re- 
spectively. For n c = 9, the Neel configuration has only a 
20% probability; the remaining 80% represent quantum 
fluctuations [see Figs. E^b) and (c)]. 

In spite of their complicated representation in the 
single-spin product basis, |0) and |1) are in many re- 
spects very similar to the states | f) and | J.) of a single 
spin and, hence, define the computational basis for uni- 
versal quantum computing with spin cluster qubits. The 



reason for this is that {|0), |1)} belong to one S = 1/2 
doublet such that 



S-|o> = |i>, S+\1) = \0), 



(9) 



where S = S x ± iS y are the spin ladder operators of the 
spin chain. 



2. Initialization 

Initialization of the spin cluster qubit can be achieved 
by cooling in a magnetic field with g^sB z <C A to tem- 
peratures T <C gfJ-BBz/ks- The spin cluster will relax 
to |0) within the spin relaxation time. Because ther- 
malization is typically slow, the preparation of a given 
initial state could be facilitated by the measurement of 
the spins within the clusters, possibly followed by local 
operations. In this way, a state which is energetically 
close to the ground state could be prepared. 



3. Quantum Gate operation 

The one-qubit phase shift gate U$, one-qubit rotation 
gate /Trot, and CNOT gate C/cnot constitute a univer- 
sal set of quantum gates. According to Eq. (JSJ, a mag- 
netic field that is constant over the spin cluster qubit 
acts on the states of the spin cluster qubit in the same 
way as on a single-spin qubit Constant magnetic fields 
B z (B x ) effect the one-qubit phase shift (one-qubit ro- 
tation) quantum gate without leakage. Because |0) and 
|1) are degenerate and separated from the next excited 
state by A, one-qubit quantum gates can be realized 
with high fidelity also by any spatially varying magnetic 
field for which |(1| Y%=i '/;/'/;/>'... •'•... ^ (U Iot ) and 
(0\Y,]=i9jtJ>BBj, z §j, z \0) ^ (U<f,), respectively. Here, 
gj is the electron spin g factor at site j. Such spatially 
varying fields can potentially cause leakage. However, if 
\gjHBBj\/A <C 1 and if all Bj are switched on and off 
adiabatically, i.e., on time scales long compared to h/A, 
the quantum gate fidelity is close to 100% as we will dis- 
cuss in detail in Sec. Ill Cl 

For the CNOT gate, one requires an exchange inter- 
action H„ between one or several spins of neighboring 
spin cluster qubits I and II which can be switched on 
and off, e.g., by electrical gates»i The simplest case of 
an exchange interaction between the outermost spins of 
neighboring clusters [Fig. Efe)] , H* — J*(t)s l n -s^ 1 trans- 
lates into an isotropic exchange interaction also for the 
two-qubit product basis. This remains true for any H* 
of the form 



Ha, 



Mt^-sf + J^t) J2(v)s)-s) +1 

3 = 1 



•« n s n -s n ' 



(10) 



Here, the factors Vj allow for a spatial variation of the 
intracluster exchange coupling constants during gate op- 
eration, where \vj\ < 1 and | J*| <C J is the limit relevant 
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for experiments. H„ will in general not only couple states 
within the two-qubit product basis {|00), |01), |10), |11)}, 
but will also lead to leakage. As long as J*(t) changes 
adiabatically and |«7*(i)| -C A for all times, leakage re- 
mains small. It should be noted that this adiabaticity 
condition also holds for single electrons in quantum dots 
where, however, the energy gap A is usually larger than 
for the spin cluster qubit. The action of H* can then be 
described by an effective Hamiltonian in the two-qubit 
product basis, 



H* = J, z (t)S l z S^ + ^^{S l+ S 11 - + S l -S ll+ ) 



-Jo(t)l, 



(11) 



where the roman numbers label the spin clusters, and 





= 4J.(t)| I <0|^ e ,,|0) 
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= Mt)h(o\ 

3=1 
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For the derivation of Eq. see Appendix lAl 
Because, for the isotropic chain, 



l^,J0)i = |i<i|£ J0)i|, 



(13) 



the coupling H* is isotropic also in the two-qubit product 
basis and acts on the states |0) and |1) of neighboring 
spin chains in the same way as an isotropic exchange 
interaction between two single spins. 

An explicit switching sequence for the CNOT gate 
based on single-qubit rotations and the unitary time evo- 
lution governed by the exchange interaction in Eq. 
is^ 

C/CNOT ~ e- i ' rS "/2 e i2 7 rn 1 .SV3 e i2 7 rn 2 .S I V3 [/jt(7r /2) 

xe" s lu^/2)e-" sl ^ 2 e-" s "' 2 e^ sl y/ 2 (14) 

for the general case where J* 2 ^ J*j_. Here, ni = 
(1, -1, 1)/V3 and n 2 = (1, 1, -1)/V3, and we have de- 
fined the unitary time evolution operator ?7*(7r/2) = 

T t exp (-i JdtH*/hj, with - / dt J*±(t)/h = tt/2. 

The gate operation time for U*(tt/2) is limited from 
below by ft/16J*|i(l|s^ iX |0)i||ii(l|s^ x |0)ii|, where J* is 
the maximum value of the exchange coupling. Matrix 
elements such as |i(l|s^ x |0)i| decrease with increasing 
n c , which leads to an increase in gate operation time. 
For realistic parameters and small n c (see Sec. Ill Fl be- 
low), J» is limited by experimental constraints rather 
than the condition J„ <C A. Then, the increase in 
gate operation time compared to single-spin qubits is 
1/(4|i(1|4 ci J0)i||ii(1|s^J0)ii|) and depends only on the 
matrix elements of the spin operators. Similarly, for a 



given magnetic field, the increase of the time required 
for a one-qubit rotation depends only on the matrix ele- 
ments |(1|sj jX |0)|. 



4- Decoherence 

For spin clusters, decoherence usually is faster than 
for single spins. The scaling of the decoherence time 
with system size depends on the microscopic decoher- 
ence mechanism. For electron spins in quantum dots, 
fluctuating fields and nuclear spins have been identified 
as dominant sourcesA^iSSiSSiSL 2 ^ In the following, we 
discuss two scenarios for decoherence in which (i) a fluc- 
tuating field couples to the total magnetic moment of the 
spin cluster qubit, and (ii) independent fluctuating fields 
act on each spin of the spin cluster qubit individually. 
In order to obtain analytical estimates for the scaling of 
with n c , we restrict our analysis to a phenomenological 
model in which decoherence is effected by a fluctuating 
classical field hb(t). In case (i), 



H B = hb(t)S z . 

The decoherence rate is obtained fromSS 
1 



(15) 



T 4> 



it ]T \(k\S z \0)\ 2 C(E k -E ) 

\k)jt\0) 

+7T £ \(k\S z \l)\ 2 C{E k -E 1 ) 
+7 r((0|5 2 |0)-(l|5 2 |l)) 2 C(0), 



(16) 



where 



1 f°° 

C{E) = — dte lEt/h (b(t)b{0)) (17) 

is the spectral density of the random field in Eq. JTSJ. 
Because of Eq. ||8J, we find that the decoherence rate 
I/t^ = 7rC(0) is independent of n c . This result, which 
is in stark contrast to the standard result that the de- 
coherence rate increases with system size, can be traced 
back to the fact that in Eq. I|15|) only the total mag- 
netic moment couples to the fluctuating field and the 
magnetic moment ±(?/is/2 of the spin cluster qubit is 
identical to the one of a single spin. Similarly, for a non- 
diagonal coupling H B — hb(t)S x , from Eq. (|16J) we find 
1/t* = 7rC(0)/2. 

Decoherence due to the coupling to nuclear spins 
is a complicated theoretical problem in its own 
right. 2 ^2^2Si2L2& In order to obtain a heuristic estimate 
for the scaling of the decoherence time with n c for this de- 
coherence scenario, we consider fluctuating classical fields 
which act independently on each site of the spin cluster, 



H 



(18) 



3=1 
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where (6,(t)6j(0)) oc §ij. For Gaussian white noise with 



(b i (t)b j (0))=2^ B 6 ij S(t), 



(19) 



I/t^ ~ n c n^ B scales linearly with n c , i.e., if fluctuating 
fields act independently on the individual spins of the 
cluster, the decoherence rate increases with n c . This re- 
sult for the scaling of l/r^ with system size agrees with 
the standard result for n c qubits which are subject to in- 
dependently fluctuating fields. 30 Note that an increase in 
the decoherence rate by a factor n c implies that the prob- 
ability for the state of the spin cluster qubit to remain 
unaffected by the fluctuating fields for a certain time t 
decreases exponentially with n c fRef. l30l) . 

Magnetic dipolar interactions are another source of de- 
coherence. Consider two single-spin qubits S\ and s 2 with 
coordinate vectors ri and r 2 , respectively, that are cou- 
pled by the magnetic dipolar interaction, 



-Hdip = 



Mo(5Mb) 2 si • s 2 - 3(si • e r )(s 2 • e r ) 



An 



d 3 



(20) 



where r = ri — r 2 , d = |r| is the distance between the 
single-spin qubits, and e r = r/|r|. In first order pertur- 
bation theory, the dipolar interaction leads to an energy 
shift ±i?dip for the product states with parallel and an- 
tiparallel spin configurations, respectively. For lateral 
quantum dots with spin quantization axis perpendicular 
to e r , 



£ di p = <T,T |#dip| T,T) 

Mg^Es) 2 (Tit |si -s 2 | t,T) _ Mg^s) 2 



An 



I6nd 3 



(21) 



because the expectation value of (si • e,.)(s 2 • e r ) vanishes 
identically. Similarly, we find 



(4, 4 l-Sdipl i, 4) — +Edip, 

(4,T l^dipl 4,T) = -E dip , 
(T,4 l^dipl T,4) = -E dip . 



(22) 



The deterministic phase shift due to the different dipolar 
energies for parallel and antiparallel spin configurations 
could be accounted for at the end of a quantum algorithm 
if no other error sources were present. However, it leads 
to correlated errors when spins are coupled to an envi- 
ronment that induces, e.g., spin flip errors. While -Edip 
is too small to induce errors of order 10~ 4 for electron 
spins in quantum dots, dipolar interactions are large for, 
e.g., P dopants in a Si matrix (see Sec. IIIFjl . We show 
next that dipolar interactions for spin cluster qubits are 
smaller than for single-spin qubits. For definiteness, we 
consider two clusters with spins at sites rj = jde y and 
rj 1 = (j + n c )de y for clusters I and II, respectively, with 
j = 1,2, ...,n c [Fig. HJc)]. The intracluster dipolar in- 
teractions lead to an unimportant energy shift that is 
identical for both |0) and |1), 



i>j 
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3s ^y s m\i) (23) 



because of symmetry. The interqubit dipolar interaction, 



7JSC 

-"dip — 



vo(gvB) 2 

And 3 
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(i - j + n c ) 3 



gives rise to an energy shift of the states 1 00) and |11) 
relative to |01) and |10), where |0) and |1) denote the 
logical basis of the spin cluster qubit, 



<00|JS£ p |00> 

(oi|^ P |oi) 



(iWplH) 



E 



dip' 



(101^1 10) = (25) 



similarly to single-spin qubits. Evaluating the matrix el- 
ements numerically, we find that the characteristic dipo- 
lar energy EJ? decreases with increasing n c We find 
E di P / E <Hp = 0.42,0.25,0.16, and 0.12 for n c = 3,5,7, 
and 9, respectively, where Edi p is the dipolar energy for 
single-spin qubits. This shows that decoherence effected 
by magnetic dipolar interactions is indeed smaller for spin 
cluster qubits than for single spins. 



5. Readout 

Readout of the spin cluster qubit is achieved by mea- 
suring the individual spins within the cluster (sj iZ ) or 
the state of the total spin of the cluster (S z ). The mea- 
surement of individual spins still is a considerable ex- 
perimental challenge. However, as has been shown the- 
oretically) 4 *^ measurement of single spins is feasible via 
charge degrees of freedom. More specifically, the state of 
a single spin on a quantum dot can be detected by a cur- 
rent flowing between spin polarized leads that are tunnel 
coupled to the quantum dot^i 

If an experimental technique is established that al- 
lows one to measure the state of a single spin, it will 
also be possible to measure the state of a spin cluster 
qubit by measurement of all spins of the cluster. Be- 
cause (0| ££1 s M |0) = 1/2 and (1| ^=1 = -1/2, 
this will allow one to unambiguously determine the state 
of the cluster qubit. However, the state of the cluster 
determines the local spin density at each site [Fig.^Jd)], 
and a probabilistic readout is possible also by measure- 
ment of single spins only. For example, for n c = 9, if the 
measurement of the central spin of the chain yields +1/2, 
the spin cluster qubit has been in state |0) with a proba- 
bility of 70%. A selective readout of several spins of the 
spin cluster qubit would also reduce the requirements on 
the readout sensitivity. For example, the sublattice spin 
(0|si, z + S3, z + S5, 2 + §r jZ + sg,z|0) ~ 1 for n c — 9 is twice 
as large as the total magnetic moment (OlS^O) and could 
be measured after separating the electron spins at even 
and odd sites of the cluster. This illustrates that readout 
of single spins is sufficient but not necessary to measure 
the state of the spin cluster qubit. 
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Spin readout can also be effected by optical means. A 
scheme for gate operation based on Pauli blocking has 
been suggested recently^ In a similar way, the creation 
of a spin polarized exciton by a circularly polarized laser 
beam and the subsequent obervation of the reemitted 
photon would allow one to determine the spin state of 
a single electron. The underlying principle is that, for a 
quantum dot with an electron spin pointing up, Pauli's 
exclusion principle prohibits the creation of an exciton 
in which the electron with spin up occupies the same 
orbital state. The creation of an exciton, which can be 
observed by the photon emitted on recombination, is pos- 
sible only if the electron spin in the quantum dot points 
opposite to the one of the exciton created by the laser 
beam. However, for qubits based on spins in a single 
quantum dot, the minimum laser focus on the order of 
the wavelength puts severe constraints on an optical read- 
out scheme. These constraints are relaxed for spin cluster 
qubits where it is sufficient to have a laser spot size with 
a diamater smaller than the size of the spin cluster qubit. 



B. Varying exchange constants 

The formation of a spin cluster qubit from an odd num- 
ber of antiferromagnetically coupled spins requires little 
control over intracluster exchange constants. Although 
both the energy gap A and matrix elements such as 
(l|s,j,a:|0) depend on the spatial variation of exchange con- 
stants, the general principle of assembling several spins 
into a cluster qubit remains valid. 

In order to demonstrate the robustness of our spin clus- 
ter qubit against a variation of exchange constants, we 
return to the isotropic spin chain but now allow for vary- 
ing fj in Eq. Because the isotropic spin chain still 
exhibits a S = 1/2 ground state doublet ?2i quantum com- 
puting is possible as discussed for fj = 1 in Sec. Ill Al 
From an experimental point of view, a priori knowledge 
of the factors fj is not necessary for quantum computing. 
Rather, the relevant matrix elements such as |(l|sj )X |0)| 
can be determined experimentally. Similarly, a quan- 
tum computer could even be assembled from spin cluster 
qubits which are not identical. 

If the exchange constants can be controlled during 
sample growth, the properties of the spin cluster qubit 
can be engineered to a certain extent. For clusters 
with centrosymmetric exchange constants, the time re- 
quired to perform the square-root of swap gate U*(ir/2) 
for given J* increases with l/|(0|si !Z |0)| 2 . For n c = 9 
and fj = 1, |(0|si, z |0)| ~ 0.18 corresponding to an in- 
crease in the gate operation time for U*(tt/2) by a factor 
(0.5/0. 18) 2 ~ 7.7 compared to the single-spin qubit for 
given J*. However, by tuning the outermost couplings 

to small values f\ = f ne —i <C minj = 2 n c -2 fjt the spin 

density at the outermost sites increases and approaches 
1/3 [see dashed lines in Fig. Hid) and Appendix iBl for 
a proof of this statement]. Although the energy gap A 
also decreases and is approximated by Jfi in this limit, 



a trade-off between the increasing matrix elements and 
the decreasing energy gap would allow one to decrease 
gate operation times compared to the chain with spa- 
tially constant exchange coupling. 

C. Leakage 

We next discuss in more detail the one qubit rotation 
gate induced by a transverse magnetic field in order to 
quantify leakage. A related analysis for the CNOT gate 
has been reported in Ref. fHi 

Because 5*10) = (ft/2)|l) and S x \l) = (ft/2) |0), a mag- 
netic field that is uniform over the spin cluster acts on 
the spin cluster qubit in the same way as on a single spin 
s = 1/2. In particular, the Hamiltonian 

H' = g^ B B x {t)S x (26) 

induces a coherent rotation from |0) to |1) without leak- 
age, as implied by quantum mechanical selection rules. 
The gate operation time for a rotation by <fi is deter- 
mined by J Q dt'giiBB x (t')/h = <f) and is identical to the 
one for single spins. In contrast, the one-qubit rotation 
effected by a spatially varying magnetic field 

n c 

H' = J2gjH B B jiX (t)sj tX (27) 

will in general lead to leakage because of finite matrix ele- 
ments (i\H'\0) and (i\H'\l) ^ coupling the compu- 
tational basis to higher excited states \i) ^ |0), |1). The 
adiabatic theorem guarantees that leakage remains small 
if the Fourier transform of Bj >x vanishes for frequencies 
larger than A/ft. Even if this adiabaticity requirement 
is not met, admixing of higher excited states to {|0), |1)} 
is controlled by the parameters gj^,BBj^ x / A and remains 
small if \gj(iBBj x \ <C A for all j. 

In the following, we concentrate on n c = 5 spins. As 
shown in Fig. E^a) and Fig. 0] a magnetic field constant 
over the cluster coherently rotates |0) into |1). This is 
also evidenced by the in-phase rotation of all spins. In or- 
der to illustrate that also a spatially inhomogeneous field 
can induce the one-qubit rotation with high gate fidelity, 
we now consider a magnetic field B^ x acting only on the 
central spin (j = 3) of the cluster. The field is switched 
on instantaneously at t = 0. For t > 0, the time evolu- 
tion is then governed by the sum of the time- independent 
Hamiltonian of the spin cluster, H, and the Zeeman 
Hamiltonian H' — g^BB^,^ x sy,, x . In Fig.^Ja), for an ini- 
tial state \ip(t = 0)) = |0), we plot the projection of the 
state \tp(t)} onto the qubit basis states for g^iBB^^j J 
0. 1 (dashed line, coinciding with the solid line on this 
scale) and 0.5 (dashed-dotted line), respectively. The 
time evolution is obtained by numerical integration of the 
Schrodinger equation. For small g^sB-i^ « A~ 0.72J, 
the spatially inhomogeneous field rotates |0) into |1) with 
high fidelity. The gate fidelity | ( 1 1 /7 ro t 1 0) | 2 decreases from 
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FIG. 3: Local spin density for all sites j of a spin chain 
with n c = 5 as a function of time (4> oc t) during the one- 
qubit rotation gate, (a) A magnetic field constant over the 
cluster, H' — g/j,BB x S x , coherently rotates |0) into |1) with- 
out leakage. Here, <j> — g/j,BB x t/h. The coherent rotation 
is evidenced by the in-phase rotation of all spins, (b)-(d) 
An inhomogeneous magnetic field, H' — g/J.BB3 tX S3, x , ef- 
fects the one-qubit rotation gate with high gate fidelity if 
\g»BB 3 , x \ < A. Here, = gfi B B 3 , x t(2\(l\s3, x \0}\)/h. The 
gate fildelity decreases from 99.8% to 93.4% and 78.5% for 
increasing g/iBBg^ — 0.1J (b), 0.5J (c), and 1J (d), respec- 
tively. In the local spin density, leakage is evidenced by high- 
frequency oscillations of neighboring spins, i.e., the excitation 
of magnons. 

99.8% {gfi B B x /J = 0.1) to 93.4% for g(j. B B x = 0.5J, 
where the typical energy scale of H' becomes compara- 
ble to A. Here, U TOt = % exp (-i /* moa dt(H + H')/h\ 

with t max = h/(4\(l\s x .3\0)\gfj,BB x ) describes the time 
evolution during a 7r-rotation. 

The decrease in gate fidelity with increasing B% x can 
be understood from the local spin densities [Figs. Efb), 
(c), and (d)]. Although only the central spin is acted on 
by B 3>x , for \gnBBz^ x \ <C A all spins of the spin clus- 
ter corotate with the central spin due to the exchange 
field. The condition |<?/iB-B3, x | "C A guarantees that the 
externally induced rotation of the central spin is suffi- 
ciently slow that all spins of the chain corotate in phase. 
For \giJ,BBs >x \ of order A, the rotation of the central 
spin induced by B^ x is too fast for the remaining spins 
of the chain to follow in phase [see, e.g., Fig. |^c) for 
short times, <fr <C 1]. The spins of the chain no longer 
rotate in phase and magnons are excited [Fig. Efb)]. 
Quantum gate operation probes the spin dynamics in 
real time and, hence, may provide new insight into the 
low energy physics of spin chains. That leakage is con- 
trolled by the parameter \g/iBB3 tX \/ A can be traced back 
to the existence of a ground state doublet. 

To illustrate this point, we next contrast our results for 
a system with antiferromagnetic exchange coupling and 
a ground state doublet with a ferromagnetic chain with 
exchange coupling J < in Eq. @. For the ferromag- 




4> 

FIG. 4: Time evolution of a state \ip(t)} with 1-0(0)} = |0) 
during the one-qubit rotation induced by a spatially con- 
stant magnetic field (solid line) and a magnetic field acting 
only on the central spin of the cluster (dashed and dashed- 
dotted lines) for function of time, cf> oc t. The 

constant magnetic field effects the quantum gate with 100% 
fildelity. In order to illustrate the effect of leakage with in- 
creasing Zeeman energy for a spatially inhomogeneous field, 
H' — g/J, B B3, x S3, x , we assume that H' is switched on instan- 
taneously at t = 0. (a) For g\iBBz^ x — 0.1J (dashed line), 
the inhomogeneous magnetic field still effects the one-qubit 
rotation with 99.8% fidelity, (b) Leakage out of the computa- 
tional basis, i-(](o|v>(t))| 2 +|(i|V>(*))| 2 ), 

remains smaller than 

0. 3%. In contrast, for gfiBBs ]X — 0.5J (dashed-dotted lines), 

1. e., comparable to A ~ 0.72J, the fidelity is only 93.4% and 
leakage is of order 7%. 

netic chain, the ground state has degeneracy n c + 1. A 
computational basis could be defined also in terms of a 
subset of the ground state multiplet for J < 0, e.g., by the 
states \S Z = —n c /2) and \S Z — —n c /2 + 1). Due to the 
n c + 1-fold degeneracy of the ground state multiplet, the 
system is more prone to leakage than the antiferromag- 
netic systems considered here because quantum mechan- 
ical selection rules no longer prevent transitions from the 
computational basis to other states in the ground state 
multiplet. 

D. XY-like chains 

Our considerations have so far been focused on an 
isotropic exchange coupling which is found, e.g., for cou- 
pled quantum dots^i However, spin cluster qubits could 
not only be realized in such artificial superstructures. 
Rather, a wide variety of spin-1/2 chains in which mag- 
netic ions are coupled by an exchange interaction has 
been synthesized^ The control of single ion spins in 
a molecule or solid in which the ions are rather closely 
spaced is even more challenging than the control of a sin- 
gle spin in a quantum dot. Defining a qubit in terms of 
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the collective state of a short spin chain will make it sub- 
stantially simpler to address qubits. However, in contrast 
to coupled quantum dots, the exchange coupling in many 
spin chains found in nature is anisotropic. 33 Hence, we 
next turn to a discussion of anisotropic chains, Jj_ =/= J z . 

For odd n c , the spectrum still exhibits a ground state 
doublet of S z eigenstates with eigenvalues ±ft/2, respec- 
tively. However, S 2 is no longer a good quantum number. 
Both for J± 3> J z (XY-like systems) and J± <C J z (Ising- 
like systems), |0) and |1) can be explicitly constructed. 

We first consider the XY model with J z = in Eq. J2J . 
By the Jordan- Wigner transformation^ the XY chain is 
mapped onto a system of nonintcracting spinless fcrmions 
on a lattice with spatially varying hopping amplitudes, 



H = - 



ra c -l 

E 



where 



exp 



ITV 



E 

k=l 



Sk, 



(28) 



(29) 



annihilates a Jordan- Wigner fermion at site j. The 
problem is thus reduced to calculating the one-particle 
eigenenergies and eigenstates of Eq. (|28p. The one- 
particle Hamiltonian has (n c — l)/2 pairs of states with 
negative and positive energy =fEi, respectively, which are 
related to each other by staggering of the wave function, 



= (ei,i,±e ij2 , ej, 3 , ±e iiS , . . .) , 



(30) 



where e^j are real numbers. In addition, there is one 



eigenstate 
e cx (l,0,~, 0,^,0, 



/l/3 ■ ■ ■ fn c -2 
fifi ■ ■ ■ fn c -l 



(31) 



with energy eigenvalue 0. The ground state doublet of the 
XY chain corresponds to the lowest (n c — l)/2 and (n c + 
l)/2 Jordan- Wigner fermion levels filled. Similarly to 
the spin chain with isotropic exchange interactions, one- 
qubit gates can be realized by magnetic fields B z (t) and 
B x (t). By numerical exact diagonalization of small spin 
chains (n c — 9), we have shown that | ^1 j^Sv^jO) | remains 
of order 1/2 for various set of fj [e.g., fj = 1, fj — 
sin(j7r/n c )], such that the operation time for the one- 
qubit rotation gate is only limited by ft/A. An isotropic 
interqubit coupling = J*(t)sJ, •§]* for two-qubit gates 
still translates into the effective Hamiltonian in Eq. 
With Eqs. I|30l) and (|31[) . from the completeness relation 



Ei:r 13/2 «? + 



'0,j 



1 for j = 1, 



and 



fyi/jj — 1/2, one can calculate all matrix elements 



°J,* - Y 3 

entering the effective coupling Hamiltonian Eq. 



c 0,l 

eo,i 



(32) 
(33) 



where eo.i is the first component of the normed one- 
particle eigenstate defined in Eq. (|31|l . In particular, 



for h 



J0|s„ c , z |0) = 1/K + l) and |(l|s nc , x |0)| - 
1). Due to the anisotropy of the intra- 



1/V2K 

chain exchange coupling, iJ» (which is isotropic in the 
single-spin operators) translates into an effective XXZ- 
Hamiltonian in the two-qubit product basis. Neverthe- 
less, the CNOT gate can still be realized according to 
Eq. 114fl . For the anisotropic chain, a magnetic field ap- 
plied along an axis n translates into a rotation around the 
axis (n x ,n y ,n z /2\(l\S x \0)\) in the Hilbert space spanned 

by{|0),|i)}, 

H' = g^BBn ■ S 

~ ffM Bi?2|(l|^|0)|[n :c (|0)(l| + |l)(0|)/2 

+ % i(|l)(0|-|0){l|)/2] (34) 
+ ff MflBn,(|0)<0|-|l)<l|)/2. 

A one-qubit rotation around an arbitrary axis, e.g. ni 
in Eq. I|14ll . hence requires appropriate rescaling of the 
magnetic fields i& 

E. Ising-like chains 

In the Ising limit J z 3> Jj_ the ground state doublet 

|0) = I t)l| 1)2... I Vn c + 0(J ± /J z ), 

|1> = I Oil T>2...|i)„ o + 0(J±/J*) (35) 

is separated from the next excited state by an n c - 
independent A ~ J z min(/j). In perturbation theory in 
J±/J z , for fj = 1, the matrix elements 



KW>| 

l(l|Sn C) x|0)| 



1 /2J, 



(»c-l)/2 



4 \ J z 
2\J Z ) 



(36) 



decrease exponentially with system size because S x and 
Sn c ,x only flip one spin within the chain. Expanding the 
states |0) and |1) in powers of J±/J z it follows that finite 
matrix element of §i jX and S x between |0) and |1) occur 
only in order (n c — l)/2 in J±/J z & Even for medium 
sized chains n c > 9 and Jj_/J z < 0.2, an isotropic inter- 
qubit coupling Hamiltonian H* translates into an effec- 
tive Hamiltonian Eq. (flip of Ising form (Fig. EJ). Be- 
cause of the long gate operation times implied by Eq. I|3t)|) 
for the one-qubit rotation and, in particular, the CNOT 
gate, only quantum computing schemes which require a 
small number of such operations 3 - appear feasible. 

In Fig. [5] for a chain with n c = 9, we compare our 
analytical results (solid lines) for the matrix elements 
|(l|5a;|0)| and |(l|si ja |0)| as functions of anisotropy with 
exact diagonalization for the chain with fj = 1 (symbols) . 
Because matrix elements of order unity imply quantum 
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FIG. 5: Transition matrix elements |(l|,§a.|0)| and |(l|si,a:|0)| 
as a function of exchange anisotropy (Jj_ 7^ J z ). The matrix 
elements determine gate operation times for the one-qubit ro- 
tation and CNOT gate. Diamonds show numerical results 
obtained for n c — 9 and spatially constant exchange inter- 
actions, fj = 1, in comparison with analytical results (solid 
lines) [see Eq. 1361 1. 



gate operation times comparable to single spins, the re- 
sults in Fig. [S] show that universal quantum computing 
based on a sequence of one-qubit rotation and CNOT 
gates is feasible for a wide range of spin cluster qubits. 



F. Experimental realization 

We illustrate the advantages of a spin cluster qubit 
formed by n c = 5 spins s = 1/2 in a one-dimensional 
array of quantum dots with diameter d ~ 50 nm. For 
realistic parameters^ (J ~ 10 Kks, J* — 2.3 Kks, 
gfisB ~ 0.7 Kks) and a magnetic field which decreases 
from its maximum values B x at the central site of the 
chain to Q.2B X at sites j — 2,4, the operation time 
for one-qubit gates increases by a factor 1/2 1 (l|0.2s2.x + 
Ss,x + 0.2s4 iX |0) I ~ 2.2 compared to a single spin. The 
operation time for the square-root of SWAP gate is in- 
creased by a factor of l/(2|(0|si^|0)|) 2 ~ 4 compared to 
the single-spin qubit. However, the operation time for the 
CNOT gate as defined in Eq. (|14|> is mainly determined 
by the single-qubit operations of the sequence. Hence, for 
the minimum operation time of the CNOT gate we find 
386 ps for spin clusters instead of 165 ps for single spins. 
The decrease of decoherence time strongly depends on 
the microscopic origin of decoherence. From the heuristic 
argument in Sec. Ill Al we find that the decoherence rate 
due to globally fluctuating magnetic fields does not scale 
with n c and is equal to the one of the single-spin qubit. 
Decoherence caused by fluctuating local fields scales lin- 
early with n c . For a spin cluster qubit with n c = 5, 
we estimate that the number of quantum gates which 
can be performed during the decoherence time decreases 
by at most a factor of 10 compared to the single-spin 



qubit. However, in contrast to single-spin qubits, mag- 
netic fields and exchange constants must be controlled 
only over length scales n c d — 250 nm and 2n c d = 500 nm, 
respectively instead of d = 50 nm. This would allow one 
to control the exchange interaction between neighboring 
clusters optically^ Note that, for the small clusters an- 
alyzed here, A is so large that neither adiabaticity nor 
the requirement that the Zeeman energy and J* be small 
compared to A provides a serious restriction on quantum 
gate operation times. 

Decoherence due to dipolar interactions is unimpor- 
tant for electron spins in quantum dots with a typical 
distance d = 50 — 100 nm, where E'dip [Eq. lf2*T)l] is of 
order 10~ n Kfcs = H/0. 25s. However, for P dopants in 
a Si matrix with a typical distance d — 100 A, i?dip = 
6.2 x 10~ 7 Kfcs = H/12fxs. Recent experiments^! show 
that magnetic dipolar interactions are indeed among the 
most dominant decoherence mechanisms for P dopants 
in a Si matrix. The characteristic time scale for dephas- 
ing by dipolar interactions, 12/is, is larger than realis- 
tic gate operation times by only a factor 10 2 — l0 3 i22iSi 
However, for a spin cluster formed by n c = 9 spins with 
a dominant nearest neighbor exchange interaction, the 
intercluster dipolar interaction energy is reduced com- 
pared to the one of single spins by a factor of 8, i.e., the 
characteristic decoherence time for dipolar interactions 
increases to lOO^is. 



III. SPIN CLUSTERS IN d > 2 

So far, our considerations have been restricted to spin 
chains. The main ideas discussed above apply to a much 
larger class of antiferromagnetic systems with uncompen- 
sated sublattices. We illustrate next that quantum gates 
are feasible also if spins s = 1/2 are arranged in a two- 
or three-dimensional cluster. For definitencss, we restrict 
our attention to an isotropic exchange interaction J > 0. 



A. Bipartite lattice 

We first consider an odd number of spins arranged on 
a bipartite lattice with the number of sublattice sites dif- 
fering by 1, e.g., a rectangular lattice with L\ x L2 sites, 
where are odd. This two-dimensional lattice exhibits 
a ground state double^2i. Similarly to the spin chain, the 
computational basis {|0), |1)} can be defined in terms of 
the S^-eigenstates of the ground state doublet. Here, S 
is the operator of total spin of the two-dimensional array 
From a spin- wave ansatz for the elementary excitations, 
the energy gap A separating the ground state doublet 
from the next excited state can be estimated as 



Jtt 



min(Li, L 2 ) 



(37) 



Because the characteristic features of the ground state 
doublet carry over from the one- to the two-dimensional 
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FIG. 6: Two-dimensional spin clusters. Each dot represents 
a single-spin qubit. (a) The spin cluster qubit scheme is read- 
ily extended from spin chains to any bipartite lattice, (b) 
Spin arrays with frustrated bonds have a highly degenerate 
ground state (fourfold degeneracy for three spins), (c) If the 
frustrated bonds are part of a larger array, the high degener- 
acy is usually lifted and a ground state doublet remains. 

spin cluster qubit, quantum computing with two- 
dimensional spin arrays on bipartite lattices is possi- 
ble with the techniques discussed in Sec. Ill Al Gate 
operation times are determined by the matrix elements 
|<0js J%2; j0> j = |(l|s_,- )a ,|0)|. For the 3 x 3-lattice shown in 
Fig. Eta), from exact diagonalization we find (0|sj iZ |0) = 
0.15 for sites j in the center of the edges, (0|sj. z |0) = 0.23 
for sites j at the corners, and (0|sj )Z |Q) = 0.17 for the cen- 
tral site of the cluster. Similarly to the spin chain, the 
ground state doublet is robust against a spatial varia- 
tion in the exchange constants as long as all exchange 
constants remain antiferromagnetic. 

B. Geometrically frustrated systems 

For non-bipartite lattices, a ground-state doublet is not 
guaranteed to emerge. The simplest example is the geo- 
metrically frustrated system of three spins s = 1/2 shown 
in Fig. Eft), 

H = J(§i • S 2 + Si • S 3 + S 2 • S 3 ) 




which has a fourfold degenerate ground state with energy 
eigenvalue E = — 0.75J. The eigenstates can be chosen 
as 

|o) = (im>-iuT»/V2, 
|i> = (im>-iUT»/V2, 

1 2) = (2|TU)-UU)-|UT))/V6, (39) 



1 3) = (2|ITT)-|TTI)-|TIT))/V6. 

As demonstrated in Ref.^1 these states could still define 
a logical qubit robust against certain sources of decoher- 
ence. However, quantum gate operation would always 
require control over single exchange interactions or lo- 
cal magnetic fields and exclude quantum computing with 
control parameters which vary slowly in space. 

Geometrical frustration does, however, not in general 
rule out the existence of a ground state doublet. In the 
more generic case that geometrically frustrated bonds 
are part of a larger system [Fig. Etc)] or the exchange 
constants in Fig. Et D ) are n °t all equal, a ground state 
doublet emerges. In this case, the logical states of the 
spin cluster qubit again can be defined in terms of the 
S z eigenstates of the ground state doublet and quantum 
gate operation is possible with magnetic fields and ex- 
change constants varying slowly over the cluster. For 
systems as shown in Fig. Et c )j m which some bonds are 
frustrated, A is usually smaller than in the case of bi- 
partite lattices. For example, for Fig. Et c )j A = 0.157J 
compared to A = 0.991J for Fig. Eta). According to the 
adiabaticity requirement, the small gap limits gate oper- 
ation times more severely for the system in Fig. Etc). 



C. Experimental realization 

Because spin cluster qubits emerge also in two- 
dimensional regular spin arrays, spin cluster qubits can 
be arranged in a plane if the positions of single spins can 
be controlled as, e.g., for lateral quantum dots [Fig. Eta)]. 
For a spin cluster qubit formed by L x L = n c quantum 
dots, A cx 1/L. Decoherence due to globally fluctuating 
fields does not increase with n c , whereas independent 
local Gaussian white noise gives rise to a decoherence 
rate cx n c . Two-dimensional spin cluster qubits are 
hence particularly interesting for qubits in which deco- 
herence is induced mainly by global rather than local 
fluctuating fields. 

More importantly, a spin cluster qubit can be defined 
even for a wide range of systems in which the positions 
and exchange constants cannot be accurately controlled. 
For P atom electron spins in a Si matrix, because of rapid 
oscillations of the exchange coupling between atoms at 
large distances, placement of atoms with lattice spac- 
ing precision is required for single-spin qubits. 38 Without 
this precision, the exchange interaction at large distances 
vanishes with a probability of 50%. In contrast, for spin 
cluster qubits formed by a small number (e.g. three) of 
P dopants located close to each other, the spin defin- 
ing the logical state of the qubit is delocalized over the 
cluster. The effective exchange coupling between neigh- 
boring qubits obtained by integration of the exchange 
interaction over the clusters is finite with a high proba- 
bility. Because the intracluster exchange interaction at 
small distances varies strongly with distance^ for a ran- 
dom arrangement of three spins, the exchange couplings 
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will differ with high probability and the system is not 
frustrated. 



within a coherent state path integral formalismi 3 ^! We 
only state the main results of our calculations here. Fur- 
ther details are given in Appendix 



IV. LARGER SPINS 

So far, our considerations have been restricted to clus- 
ters formed by spins s = 1/2. We next consider antifer- 
romagnetic systems with larger spins s > 1/2. 

A. Antiferromagnetic molecular clusters 

Only very recently it has been shown theoretically that 
Grover's algorithm can be implemented with ferromag- 
netic molecular magnets using a unary encoding^i 4 ^ 
In view of universal quantum computing, ferromagnetic 
clusters such as the molecular magnets Mni2 and Fe§ 
(Ref. |43J) suffer from the large net spin which usually 
means large matrix elements coupling the spin to the en- 
vironment and, hence, short decoherence times. 

In contrast, in antiferromagnetic systems, such as var- 
ious ring compounds^ 3 " the spins couple such that they 
form a small total magnetic moment. Antiferromag- 
netic clusters which have unequal sublattice magnetiza- 
tion, will in general have a ground state multiplet rather 
than the singlet found for systems with compensated sub- 
lattice spins4^ Several antiferromagnetic molecular mag- 
nets comprised of spins with quantum numbers larger 
than 1/2 have been synthesized to dato 45 i 46 ' 47 i 48 i 49 in- 
cluding several compounds with uncompensated sublat- 
tice spins^fi As a paradigm, we consider systems with 
isotropic exchange interaction J, but allow for an easy or 
hard axis single-spin anisotropy, 



H = Jsi ■ s 2 + k z {s\ 



(40) 



Here, s± and si = Si ± 1/2 are the spin quantum num- 
bers of the two sublattices, respectively, J > is an effec- 
tive exchange constant, and k z the single ion anisotropy. 
Equation (|4(J|) has a ground state doublet {|0),|1)} of 
S z = si,z+S2,z eigenstates with eigenvalues ±fi/2, respec- 
tively. Because [S x , H] = for the Hamiltonian Eq. (f^ . 
the logical qubit basis states have an expansion of the 
form (for s 2 — Si — 1/2) 



|0) 



Si 

E 

mi- — + 1 



ft, 



ll |mi, 1/2 - mi) 



a 



- Sl +i\ - si ■+ 
" Sl |si, -si 4 



Ml -1/2) 
+ 2, Sl -3/2) + 
1/2) 



(41) 



in terms of the spin product basis. For s^i 3> 1, an- 
alytical expressions can be derived both for the action 
of a magnetic field (one-qubit rotation gate) and for the 
action of an interqubit coupling Hamiltonian (two-qubit 
gates) between clusters I and II, 



B. Hard axis systems 

For strong hard axis anisotropy, k z > and 
k z s\ 2 / J 3> 1, the spins Si and S2 lie close to the x-y 
plane for both states of the ground state doublet. A large 
contribution in the expansion Eq. I|41|) comes from the 
states | mi = 0)|m 2 = 1/2) and |mi = 0)|m 2 = —1/2), 
respectively. For illustration, for s\ — 3, S2 = 5/2, and 
k z /J = 0.2, by numerical diagonalization of Eq. i|4TJ|) . we 
find 



10) = 



0.25 
+0.52 
+0.42 



3,--) -0.41 



2,- 



!,-->- 0.52 



(43) 



■1, 



0.24 



The state |1) is obtained by \m\, mi) — » | — mi, —m^) on 
the right hand side of Eq. I|43|l . In agreement with the 
semiclassical theory, a major contribution to {|0),|1)} 
comes from states with small mi and m%. 

In the following, we restrict our attention to systems 



with large anisotropy, k z (s\ - 
(Appendix EJ) and 

KW)| 

l(l|*i,»|0)| 

l(i|^|o)| 



f)/J > 1. Then, A ~ J 



1/4, 

si/2, 

sill. 



(44) 



In particular, Eq. (|42|) translates into the effective Hamil- 
tonian 



= J*|(0|s M |0)||(0|s 2 ,,|0)| 



+ - 



J*SiS2 



/I 




Vo 

/ 
10 
10 

Vo o o o 






-1 





(45) 



(42) 



in the two-qubit product basis. As discussed in Sec. Ill Al 
the CNOT gate can be realized with a unitary time evo- 
lution governed by this effective qubit coupling of the 
XXZ-form. 

Matrix elements of order unity in Eq. (|44|l show that, 
e.g., a magnetic field B x efficiently rotates the state |0) 
into |1). This is not a priori evident given the rather 
complicated representation of the ground state doublet 
in the single-spin product basis [Eq. (|4ip]. The large 
matrix elements arise because, for both |0) and |1), the 
spins lie close to the x-y plane in the hard axis system. 
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the electrical conductivity of individual molecules can be 
switched between two states in a controlled way. By 
connecting molecular magnetic clusters by reversible re- 
dox switches one could also switch intercluster exchange 
paths. Alternatively, if the relative position of molecu- 
lar magnetic clusters can be controlled, the intercluster 
exchange interaction can be switched on and off via the 
overlap of electron orbital wave functions by moving clus- 
ters relative to each other. 



CONCLUSION 



FIG. 7: Matrix elements of spin operators of a qubit formed 
by two spins with spin quantum numbers Si = 7 and S2 = 6.5. 
Numerical data (symbols) obtained from exact diagonaliza- 
tion are in good agreement with analytical results (solid lines). 

C. Easy axis systems 

For k z < 0, configurations with spins aligned along the 
z-axis are energetically favorable. We restrict our atten- 
tion to systems with large anisotropy, A\k z \{s\+ s\) / J ^> 
1. Because a transition from |0) to |1) requires a rotation 
of both spins through a large energy barrier, from the the- 
ory of spin quantum tunneling in antiferromagnetic sys- 
temaa Swc find that Kl |&lO)|, |(l|si |3; |0)|, |(l|a 2 ,x|0)| oc 
exp(— \J&\k z \(s\ + s^jjj) <C 1 are exponentially small. 
Similarly to a spin chain in the fsing limit [Sec. Ill E] . the 
easy axis system is a candidate for quantum computing 
schemes as suggested in Ref. |3^ 

The analytical results for the matrix elements dis- 
cussed here are compared with numerical exact diago- 
nalization for si = 7 in Fig. [7| We find good agreement 
with our semiclassical results. 



D. Experimental realization 

Single molecule electrical switches^^^S have nour- 
ished hopes that, in the future, it will be possible to 
down-scale computers to the level at which bits or qubits 
are represented by single molecules. The results in 
Sees. IIVBI and IIV CI show that, in such bottom-up ap- 
proaches aiming at a universal quantum computer, con- 
trol is not required at the level of single atom spins but 
only on the scale of molecule spins. In particular, molec- 
ular magnetic clusters with an effective spin S — 1/2 de- 
fine a qubit. One-qubit quantum gates could be effected, 
for example, by a magnetic tip as used in magnetic force 
microscopy^ The spatial resolution of these techniques 
currently lies in the range of 10 — 100 nm (Ref. l59|) and 
approaches the typical size of molecular magnetic clus- 
ters^ 

Control of the exchange interaction between molecules 
is challenging. As has been demonstrated recently^^ 



In conclusion, we have shown that quantum computing 
is possible with a wide variety of clusters assembled from 
antiferromagnetically coupled spins which form an effec- 
tive total spin S = 1/2. For arrays of spins s = 1/2, the 
existence of a spin cluster qubit requires little control over 
the placement and intracluster coupling of the spins and 
the spatial dimension of the array. This remains true for 
a wide range of systems with uncompensated sublattice 
spins differing by 1/2. We have shown explicitly that, 
for the spin cluster qubit, initialization, quantum gate 
operation, and readout are possible with the techniques 
proposed and analyzed for single spins. The scaling of 
the decoherence time with system size strongly depends 
on the microscopic decoherence mechanism. Spin cluster 
qubits are particularly promising in situations where de- 
coherence is induced mainly by globally fluctuating fields 
during quantum gate operation and the decoherence rate 
of the spin cluster qubit is comparable to that of a single- 
spin qubit or for systems in which magnetic dipolar inter- 
actions are the dominant decoherence mechanism. The 
main advantage of spin clusters compared to single spins 
is that requirements on local control of magnetic fields 
and exchange interactions can be traded for longer gate 
operation times. We have illustrated the feasibility and 
advantages of spin cluster qubits for arrays of quantum 
dots, P dopants in a Si matrix, and molecular magnetic 
clusters. 

In contrast to single spins, spin clusters are not intrinsi- 
cally two-state systems and leakage during quantum gate 
operation must be accounted for. For the one-qubit ro- 
tation gate, we have shown that leakage is small if the 
magnetic field which induces the rotation is switched on 
and off adiabatically or if the Zeeman energy remains 
small compared to A. 

Finally, we note that, because any qubit can be 
mapped onto a spin s = 1/2, the results of this work do 
not only apply to quantum computing proposals based on 
spin degrees of freedom but to any quantum computing 
scheme. More specifically, for any qubit for which meth- 
ods for initialization, quantum gate operation, quantum 
error correction, and readout have been identified, a clus- 
ter qubit can be formed by coupling several qubits. For 
the cluster qubit, initialization, quantum gate operation, 
quantum error correction, and readout are possible using 
the same techniques as for the original qubit. 
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APPENDIX A: EFFECTIVE EXCHANGE 
HAMILTONIAN 

Here, we derive explicitly the effective coupling Hamil- 
tonian in Eq. (|ll|l from Eq. (|10|l . The first and second 
term in Eq. I|ll|) result from J*(i)sJ^ ■ s^ 1 in the micro- 
scopic coupling. Decomposing • sP = s l n z sY z + 



n c -l 



(4t 



I+£lI- 



si-sii 



s 1 + )/2 in terms of spin ladder opera- 



tors, one can readily evaluate the matrix elements in 
the two-qubit product basis. Because, by definition, 
S Z \Q) = (ft/2)|0) and S z \l) = (-ft/2)|l) and a^.S*. 
conserves the z-component of total spin in each cluster 
separately it follows that 



a si 



11)110)! 



0. 



(Al) 



Similarly, all other off-diagonal elements of J>\ z van- 
ish. Because of 



£|o>!= !(i|4tli)i = o 



(A2) 



i+s"- 



the transverse exchange interaction J*(i)(s„+s 
s^~s^ I+ )/2 has finite matrix elements only between the 
states |0)i|l)n and 1 1) 1 1 0)n . This completes the proof 
that the intercluster exchange term J*(f)sjj ■ s^ 1 leads to 
the first and second term in Eq. i|ll|). 

It remains to show that a possible change in intraclus- 
ter exchange interaction constants during two-qubit gate 
operation, J* (t) J2]=i (Vjit) ■ s) +1 + vfsf ■ sf +1 ), only 
leads to a term proportional to 1 in Eq. i|ll|) . This term 
conserves all components of the total spin of clusters I 
and II, 



[Si £ v]s]-s] +1 } = 0, 

3=1 



(A3) 



for a — x,y, z, and similarly for II. Hence, all off-diagonal 
matrix elements such as i(l|ujsj- ■Sj+i |0)i vanish. Finally, 

because |l)i = S 1 '^! = 2£*|0)i, with Eq. l|A3j) . 



71c — 1 



x (l| n<0| J.(t) J2 ^•^•+i|l>i|0>n 



3 = 1 



Tic — 1 



i<0| n <0|25i J,(t) J2 v)s) ■ s^SllOmM) 

3 = 1 



= i<0| n <0| J.(t) J2 v }s) •§J + i(2^) 2 |0)i|0) II 



3=1 
n c -l 



i(0| „(0| ,h(t) ^sj.. s i +1 |0)i|0)n. 

3 = 1 



In the second line of Eq. (|A4|) we have invoked that 
{|0)i, |l)i} belong to one spin-1/2 doublet, the third line 
then follows from Eq. I|A3[) . With a similar argument it 
can be shown that all diagonal matrix elements in the 
two-qubit product basis are equal and J*(t) $27=1 i v ^} ' 
Sj +1 + uPsP • s? +1 ) translates into a term J (t)l in the 
effective coupling Hamiltonian Eq. I|ll|) . 

Finally, we prove Eq. (|13fl which implies that, for 
isotropic intracluster exchange interactions, the effective 
two-qubit Hamiltonian is also of Heisenberg form. For 
simplicity, we omit the label I of the spin cluster qubit in 
the following. In order to formally calculate (0|s„ ciZ |0) 
and |(l|s nt , jX |0)|, we define the spin operators 



S~ — S 



n c — l 

E 

3 = 1 



'J, a 



(A5) 



of all but the outermost spin j = n t 
erally, |0) can be expanded as 



of the cluster. Gen- 



|0) = a |*)|t) + fe|$)|4), 



(A6) 



where \i$>) and |$) describe the normed states of the left- 
most n c — 1 spins in the array and a and b are real 
numbers. Because S z \0) = (K/2)\0), |*) and |$) are 
eigenstates of S' z with eigenvalues and ft, respectively. 
S x \0) = (ft/2) 1 1) is an S z eigenstate with eigenvalue 
-(ft/2), such that 6|$) = -aS' + \^), and 

|0) = a(|*)|T)-^|*)U», 

|1> = a(^_|*)|T>-(l-^_i;)|*>U», 

1 (A7) 



'1 + (*|5LS , ;|4') 
From Eq. (|A7fl we calculate 

(0|s„ CiZ |0) = -{l\s neiZ \l) = -- 



2i + ms'_s',m 



(l|Sn c ,x|0) 



1 1 - ms'_ s\ m 



2l + (V\S'_ S' + \^} 
which proves Eq. I|13|) . 



(A8) 



APPENDIX B: ISOTROPIC CHAIN WITH 
SPATIALLY VARYING EXCHANGE 
INTERACTION 

The local spin density in the energy eigenstates of 
Eq. depends sensitively on spatial variations of the 
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exchange interaction [Fig. |2{d)]. Whereas for fj = 1 
(solid lines) the magnetization density in each of the 
sublattices increases toward the center of the chain jEtt 
the opposite behavior is observed for an exchange in- 
teraction fj = sin(j7r/n c ) (dashed lines). In the limit 
fx = fn c -i *C minj = 2 ) ..., rae _2 fj, the increase of local 
spin density toward the ends of the chain can be under- 
stood quantitatively. The ground state doublet of the 
spin cluster qubit can be constructed explicitly from the 
ground state doublet {|0)„ c _2, |l)n c -2} 01 the chain with 
the outermost spins removed. For Jfi much smaller than 
the energy gap A„ c _2 of the chain formed by the n c — 2 
central spins, the coupling of the outermost spins can be 
treated perturbatively. For the chain with centrosym- 
metric exchange couplings, fj = fn c -j, from the ansatz 



|0) = ai| T)|l)„ c - 2 | T)+a 2 | T)|0)„ c - 2 | i) 
+a 3 \ i)|0) We - 2 | T) + 0(J/i/A nc _ 2 ), 



(Bl) 



we find (ai, a 2 , a^) = 
state of the chain with n 

lim 



'2,— 1,— 1)/V6 for the ground 
; spins and, hence, 



|(0|s 1>z |0)| = l/3. 



(B2) 



APPENDIX C: LARGE SPINS 



-iAshd>(l — cost 



(C2) 



where the last factor accounts for the difference As = 
32 — si of the spin quantum numbers and <j) — d T (f> is 
the imaginary time derivative. In the limit of strong 
anisotropy, k z (sf + 82)/ J ^> 1, Eq. I|C2|I can be expanded 
to second order in 9 — tt/2 and the fluctuations are inte- 
grated out, leading to 



Z = / T>4> exp - 



h8 



dTL E [<p}/h 



(C3) 



with an effective Euclidean Lagrangean 



^2 ±2 



2J 



iAs he 



(C4) 



After continuation to real time, by a canonical transfor- 
mation we obtain the Hamiltonian of the system in terms 
of the Neel vector operator, 



(C5) 



The matrix elements in Eq. (|44l) can be calculated from 
coherent state spin path integrals?^ We focus on strong 
easy plane systems, k z > and fc z (sf + s§)/ J ^> 1. Fol- 
lowing the standard approach for antiferromagnetic sys- 
tems, the partition function of the two-spin system is 
expressed as path integral over the Neel vector n and ho- 
mogeneous magnetization 1 defined by Si = sin + 1 and 
S2 = — s 2 n + 1, where n • 1 = 0. Integrating out 1 in a 
saddle point approximation and parametrizing 



(CI) 



the Euclidean action of the system can be written as^SS 



where 




L E = — 
2J 



+ sin^ 9 <jf) + fc z (sf + S2 ) cos^ 



COS( 

sin < 




, 1 = 






pf/h - As 



(C6) 



and pij) is the momentum operator conjugate to the in- 
plane polar angle <f>, \p$, <f>] — —ih. By inspection of 
Eq. i|C5(l . we find that the spin system Eq. (|40|) has 
been mapped onto the Hamiltonian of a particle on a 
ring threaded by a magnetic flux cx As. In particu- 
lar, for half-integer As, the Hamiltonian has a ground 
state doublet {|0), |1)} with wave functions ipo(<l>) — 
exp(i(m+ l)<j))/*\/2~n and ipi((f>) = exp(im(j>) / ^/2tt , where 

m = [As\. From §i jX — sicos^ and §2,x — — S2COS0, 
one immediately obtains Eq. 1 14(1 . 
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